r  AGENCY  0S«  ONLY  (Uai*  bUnk)  I  2^ROAT^OATj  ^ 


4.  TTTU  AND  SUBTITLE 

ADVANCED  TECHNIQUES  FOR  ARRAY  PROCESSING 
FINAL  REPORT 


nrTniTr-FOi 


Benjamin  Friedlander 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS<ES) 

SIGNAL  PROCESSING  TECHNOLOGY,  LT! 
703  coastland  Drive 
Palo  Alto,  CA  94303 


DAALOl-fl-L-OOO  7 


t&K.  -  ' 

laiL 

.  %  J’JLi  : 

at 

y 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  ANO  ADORESS<ES) 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


10.  SPONSORING  /MONITORING 
AGENCY  REPORT  NUMBER 


/\£o 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author (a)  and  should  not  be  construed  as  an  official  Departaent  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  docunentatlon. 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  12k  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


IS.  ABSTRACT  (kktjamum 

Array  processing  technology  is  expected  to  be  a  key  element  in  communication  systems  de¬ 
signed  for  the  crowded  and  hostile  environment  of  the  future  battlefield.  While  advanced 
array  processing  techniques  have  been  under  development  for  some  time,  their  practical  nse 
has  been  very  limited.  This  project  addressed  some  of  the  issues  which  need  to  be  resolved 
for  a  successful  transition  of  these  promising  techniques  from  theory  into  practice.  The 
main  problem  which  was  studied  was  that  of  finding  the  directions  of  multiple  co-channel 
transmitters  from  measurements  collected  by  an  antenna  array.  Two  key  issues  related  to 
high-resolution  direction  finding  were  addressed:  effects  of  system  calibration  errors,  and 
effects  of  correlation  between  the  received  signals  due  to  multipath  propagation.  A  number 
of  useful  theoretical  performance  analysis  results  were  derived,  and  computationally  effi¬ 
cient  direction  estimation  algorithms  were  developed.  These  results  include:  self-calibration 
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Forward 


This  report  summarizes  the  work  performed  on  the  project  ‘‘Advanced  Techniques  for 
Array  Processing”  under  Contract  No.  DAAL03-89-C-0007,  during  the  period  from  1  March 
1989  to  30  April  1991.  The  project  considered  the  problem  of  finding  the  directions  of  multi¬ 
ple  co-channel  transmitters  from  measurements  collected  by  an  antenna  array.  We  addressed 
two  key  issues  related  to  high-resolution  direction  finding:  effects  of  system  calibration  er¬ 
rors,  and  effects  of  correlation  between  the  rece’’  ed  signals  due  to  multipath  propagation. 
This  final  report  provides  a  brief  overview  of  the  results,  deferring  the  details  to  previous 
project  publications. 
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1.  Statement  of  the  Problem 


Reliable  communications  in  the  crowded  and  hostile  environment  of  the  future  battlefield 
is  a  difficult  and  challenging  problem.  Array  processing  technology  must  be  a  key  element 
in  any  system  designed  to  face  these  challenges.  While  some  of  the  technologies  needed 
for  operating  in  this  difficult  environment,  such  as  spread  spectrum  communications  and 
advanced  modulation  and  coding  techniques,  are  fairly  mature  and  well  developed,  array 
processing  is  a  relatively  new  technology  and  its  potential  for  providing  significant  gains  in 
system  performance  has  not  been  fully  realized  yet.  Research  and  development  efforts  in  the 
area  of  array  processing  are,  therefore,  expected  to  have  relatively  high  payoff  during  the 
next  decade. 

While  advanced  array  processing  techniques  have  been  discussed  in  the  literature  during 
the  past  ten  years,  the  practical  implementation  of  these  techniques  is  still  in  an  early  stage 
of  development.  Various  experimental  systems  using  techniques  such  as  MUSIC  [41]  are  now 
in  existence,  but  as  far  as  we  know,  at  this  time  there  are  no  fielded  systems  employing 
high-resolution  array  processing  algorithms. 

In  this  project  we  addressed  some  of  the  many  open  issues  which  need  to  be  resolved 
for  a  successful  transition  of  advanced  array  processing  techniques  from  theory  into  practice. 
In  particular,  we  considered  the  problem  of  finding  the  directions  of  multiple  co-channel 
transmitters  from  measurements  collected  by  an  antenna  array.  We  addressed  two  key  issues 
related  to  high- resolution  direction  finding:  effects  of  system  calibration  errors,  and  effects 
of  correlation  between  the  received  signals  due  to  multipath  propagation. 


2.  Summary  of  the  Key  Results 

The  results  of  this  work  are  summarized  in  40  project  publications  which  are  listed  in 
section  3.  Here  we  present  only  a  brief  overview  of  this  work.  References  [1]— [40]  correspond 
to  the  project  publications  listed  in  section  3,  while  references  [41]— [48]  are  listed  in  the 
bibliography  section  at  the  end  of  the  report. 


2.1  Calibration  Issues 

A  focal  point  for  our  research  has  been  the  question:  how  to  apply  advanced  high- 
resolution  array  processing  techniques  in  the  presence  of  various  uncertainties  in  the  array 
parameters.  These  uncertainties  can  take  on  many  different  forms:  imprecisely  known  loca¬ 
tion  of  the  airay  elements,  unknown  variations  in  the  gains  and  phases  of  individual  elements, 
and  other  deviations  between  the  model  on  which  these  techniques  are  based,  and  reality. 

The  problem  outlined  above  is  of  considerable  practical  importance.  The  accurate  cali¬ 
bration  of  the  array,  required  in  order  to  apply  “standard”  high- resolution  techniques,  is  very 
costly.  Furthermore,  in  many  situations  it  is  practically  impossible  to  maintain  the  array 
calibration.  Consider  for  example  a  mobile  communication  antenna  array  which  needs  to  be 
dismantled  and  reassembled  in  the  field,  its  calibration  changing  each  time  it  is  moved.  Or 
consider  a  conformal  array  mounted  on  the  wing  of  an  aircraft,  where  mechanical  vibrations 
constantly  change  the  shape  of  the  array. 

To  address  these  problems,  which  are  a  major  obstacle  on  the  road  to  a  practical  applica¬ 
tion  of  these  promising  high- resolution  techniques,  we  have  developed  a  number  of  techniques 
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to  self-calibrate  the  array  using  “signals-of-opportunity” .  We  were  able  to  show  that  it  is 
possible  to  simultaneously  estimate  unknown  source  locations  and  the  unknown  array  pa¬ 
rameters,  under  some  reasonable  assumptions  (3,  6,  10,  11,  18,  39,  26].  It  is  important  to 
note  that  this  self-calibration  procedure  does  not  require  that  we  know  the  directions  of  the 
calibrating  sources. 

More  recently,  we  performed  an  extensive  sensitivity  analysis  of  direction  finding  tech¬ 
niques  based  on  the  covariance  matrix  of  the  received  signals  [12, 13,  25].  Using  the  results 
of  this  analysis  we  can  now  quantify  the  effect  of  various  types  of  calibration  errors  on  the 
accuracy  of  the  direction-of-arrival  (DOA)  estimates,  and  on  the  ability  of  the  array  to  re¬ 
solve  closely  spaced  sources.  These  results  appear  to  be  very  useful  for  the  proper  design  of 
high-performance  direction  finding  systems.  These  results  also  indicate  clearly  the  extreme 
ill-conditioning  of  the  problem  of  resolving  multiple  co-channel  signals  which  are  oeperated 
by  less  than  a  beamwidth  [19]. 

Our  work  on  the  sensitivity  issue  and  the  self-calibration  problem  appears  to  be  the 
first  systematic  attempt  to  analyze  and  improve  the  robustness  of  high-resolution  direction 
finding  techniques  in  the  presence  of  various  system  errors. 


2.2  Direction  Finding  for  Correlated  Signals 

Direction  finding  in  the  presence  of  coherent  multipath  is  a  difficult  problem.  The  mul¬ 
tipath  causes  an  apparent  distortion  of  the  wavefront  of  the  propagating  signal  which  can 
introduce  large  errors  into  the  DOA  estimates  produced  by  conventional  direction  finding 
systems.  High-resolution  direction  finding  systems  attempt  to  alleviate  this  problem  by  re¬ 
solving  the  components  of  the  multipath  arrivals  and  estimating  their  DOAs.  However,  there 
are  various  difficulties  associated  with  this  procedure  as  discussed  below. 

In  our  work  we  attempted  to  address  this  difficult  problem  from  several  points  of  view: 
(i)  we  studied  the  limits  on  direction  finding  accuracy  in  the  presence  of  multipath;  (ii)  we 
developed  some  practical  direction  finding  algorithms  based  on  the  eigenstructure  of  the 
covariance  matrix  of  the  vector  of  signals  received  by  the  array;  (iii)  we  investigated  the 
potential  advantages  of  using  arrays  with  polarization  diversity  to  alleviate  the  multipath 
problem.  This  work  is  summarized  briefly  in  the  following  paragraphs: 


Bounds  on  Direction  Finding  Accuracy 

In  [31],  [36]  we  Btudied  the  effects  of  correlation  on  DOA  estimation  accuracy  by  evalu¬ 
ating  the  corresponding  Cramer-Rao  bound.  We  started  by  deriving  the  Cramer- Rao  bound 
for  Gaussian  signals  with  an  unknown  covariance  matrix  (which  provides  as  a  by  product 
the  CRB  for  Gaussian  signals  with  a  knowD  covariance  matrix).  We  were  able  to  obtain 
compact,  closed  form  expressions  for  the  Fisher  Information  Matrix  associated  with  the 
joint  estimation  of  the  DOAs,  the  signal  covariance  matrix,  and  the  noise  variance. 

By  studying  the  accuracy  of  the  DOA  estimates  (as  predicted  by  the  CRB)  as  function 
of  the  magnitude  and  phase  of  the  correlation  coefficients,  we  gained  insights  into  the  effects 
of  correlation  on  direction  finding  performance.  As  expected,  DOA  accuracy  decreases  as 
the  magnitude  of  the  correlation  coefficient  increases.  We  have  also  observed  that  the  phase 
of  the  correlation  coefficient  can  have  a  strong  effect  on  DOA  accuracy,  especially  for  small 
aperture  arrays.  In  a  multipath  environment,  the  phase  of  the  correlation  coefficient  is  a 
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highly  variable,  unpredictable  quantity,  determined  by  the  difference  between  the  propaga¬ 
tion  delays  in  the  direct  and  secondary  paths.  Thus,  the  best  achievable  performance  of  a 
direction  finding  system  operating  in  a  multipath  environment  may  vary  significantly  from 
time  to  time,  due  to  fluctuations  in  the  correlation  phase. 


Direction  Finding  Algorithms 

A  number  of  direction  finding  algorithms  for  correlated  signals  have  been  proposed. 
These  include:  spatial  smoothing  algorithms  based  on  the  eigendecompoeition  of  tne  covari¬ 
ance  matrix  of  the  received  signals  [42,  43];  and  algorithms  based  on  maximum  likelihood 
estimation  [44,  45]. 

Direction  finding  techniques  based  on  the  eigen-decomposition  of  the  covariance  matrix 
of  the  vector  of  signals  received  by  an  array  of  sensors,  have  received  considerable  attention 
in  recent  years.  Their  main  advantage  over  maximum  likelihood  based  techniques,  is  their 
relative  computational  simplicity.  The  spatial  smoothing  algorithm  involves  a  one  dimen¬ 
sional  search,  compared  to  the  multi-dimensional  search  inherent  in  the  maximum  likelihood 
technique.  Using  the  root-MUSIC  approach  avoids  the  search  completely  (or  rather,  replaces 
it  with  a  root  solving  algorithm),  reducing  even  further  the  computational  complexity  of  the 
direction  finding  algorithm.  The  main  drawback  of  the  spatial  smoothing  technique  is  that 
its  use  to  restricted  to  linear  uniformly  spaced  arrays. 

In  [33,  37]  we  were  able  to  generalize  the  spatial  smoothing  technique  to  arbitrary  array 
geometries  by  using  an  interpolated  array  approach.  This  approach  is  based  on  estimating 
the  outputs  of  a  “virtual  array"  from  the  outputs  of  the  real  array.  These  outputs  are 
obtained  by  a  straightforward  linear  interpolation  technique,  with  interpolator  coefficients 
selected  so  as  to  minimize  the  interpolation  error  for  a  signal  arriving  from  a  given  sector 
(i.e.  a  range  of  bearing  angles).  Different  sets  of  interpolator  coefficients  will  be  used  to 
provide  good  estimates  for  different  sectors.  The  size  of  the  sector  will  be  chosen  to  give 
sufficiently  good  estimates  of  the  virtual  array  outputs.  The  design  of  the  interpolator  needs 
to  be  performed  only  once,  off-line. 

Given  measurements  of  the  real  *»rray  outputs,  the  interpolator  can  be  used  in  principle 
to  compute  the  outputs  of  the  virtual  array,  (more  precisely,  several  sets  of  outputs  are 
computed,  one  per  sector).  Any  direction  finding  technique  can  then  be  applied  to  the 
outputs  of  the  virtual  rather  than  the  real  array.  Actually,  it  is  not  necessary  to  compute  the 
outputs  of  the  virtual  array,  since  it  is  possible  to  compute  directly  the  sample  covariance 
matrix  corresponding  to  these  outputs.  TIk;  techniques  considered  here  require  only  the 
sample  covariance  matrix,  not  the  “raw"  data. 

If  the  interpolated  array  is  linear  and  uniformly  spaced,  the  direction  of  arrival  estimates 
can  be  computed  by  using  the  spatial  smoothing  method  in  its  root-MUSIC  version,  which 
is  computationally  simple  and  is  more  robust  to  certain  modeling  errors. 

In  [32,  34]  we  carried  out  a  performance  analysis  for  the  interpolated  spatial  smoothing 
algorithm.  Closed  form  expressions  for  the  covariance  matrix  of  the  DOA  estimation  errors 
are  derived  using  a  perturbation  analysis.  Evaluating  these  expressions  for  specific  cases  and 
comparing  them  to  the  Cramer  Rao  lower  bound  for  the  DOA  estimates,  provides  insight  into 
the  statistical  efficiency  of  this  algorithm.  The  formulas  for  the  error  covariance  are  quite 
general,  and  can  be  specialized  to  provide  results  for  other  DOA  estimation  algorithms  as 
well.  The  performance  analysis  shows  that  the  spatial  smoothing  estimator  (interpolated  or 
not)  is  not  statistically  efficient,  but  it  provides  acceptable  performance  under  a  wide  range 
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of  operating  conditions .  The  interpolated  spatial  smoothing  algorithms  were  also  tested 
by  Monte-Carlo  simulations  and  was  shown  to  perform  well  for  different  test  cases.  These 
simulations  confirmed  the  analytical  performance  predictions. 

In  summary,  the  proposed  technique  combines  the  computational  efficiency  associated 
with  eigen-structure  based  techniques  such  as  root-MUSIC,  with  good  estimation  accuracy. 
It  appears  to  be  a  very  practical  approach  to  the  problem  of  DOA  estimation  in  the  presence 
of  multipath. 


Direction  Finding  with  Diversely  Polarized  Array  s 

Antenna  arrays  with  diverse  polarization  have  some  inherent  advantages  over  uniformly 
polarized  arrays,  since  they  have  the  capability  of  separating  signals  based  on  their  polariza¬ 
tion  characteristics.  Polarization  diversity  has  been  used  quite  successfully  in  radar  systems 
to  improve  clutter  rejection  and  to  discriminate  between  different  types  of  targets.  Diversely 
polarized  arrays  are  also  useful  in  the  context  of  communications  systems. 

One  would  expect  that  polarization  diversity  would  be  quite  useful  in  high-resolution  di¬ 
rection  finding  systems  as  well.  It  is,  therefore,  somewhat  surprising  that  very  little  attention 
has  been  paid  to  this  issue  in  the  extensive  literature  on  high-resolution  array  processing. 
In  fact,  with  the  exception  of  [46],  [47]  which  present  the  extension  of  MUSIC  and  Maxi¬ 
mum  Likelihood  to  diversely  polarized  arrays,  there  seem  to  be  no  results  on  high-resolution 
direction  finding  with  these  arrays. 

Recently,  we  developed  a  computationally  efficient  algorithm  for  the  simultaneous  esti¬ 
mation  of  the  DO  As  and  polarizations  of  diversely  polarized  signals  [48].  The  algorithm  is 
based  on  polynomial  rooting  and  therefore  it  exhibits  enhanced  resolution  and  reduced  com¬ 
putational  load  (similar  to  that  of  the  root-MUSIC  algorithm  which  is  used  in  the  context 
of  uniformly  polarized  signals).  Since  rooting  methods  can  only  be  applied  to  certain  array 
configurations,  we  use  linear  interpolation  in  order  to  extend  the  technique  to  a  larger  class 
of  arrays.  Extensive  computer  simulations  were  used  to  demonstrate  the  effectiveness  of  the 
new  technique. 

In  [15,  35]  we  studied  the  performance  of  direction  finding  systems  employing  diversely 
polarized  arrays,  in  a  more  general  setting.  In  particular,  we  were  interested  in  comparing 
the  performance  of  arrays  with  and  without  polarization  diversity,  to  try  and  quantify  the 
performance  advantages  which  can  be  achieved  by  using  polarization.  Our  study  focused  on 
the  statistical  accuracy  of  the  direction  estimates,  using  the  Cramer  Rao  Bound  (CRB)  for 
the  joint  estimation  of  the  directions  of  arrival,  polarization  parameters,  signal  covariance 
matrix,  and  the  noise  variance.  We  were  able  to  develop  compact  closed  form  formulas  which 
show  the  dependence  of  the  bound  on  various  system  parameters. 

The  CRB  was  derived  under  the  assumption  that  the  signal  parameters,  namely  the 
polarization  parameters  and  signal  covariance  matrix  (as  well  as  the  variance,  of  the  mea¬ 
surement  noise)  are  unknown  and  need  to  be  estimated  as  part  of  the  direction  finding 
process.  This  bound  is  more  realistic  than  the  one  which  assumes  that  the  signal  (and  noise) 
parameters  are  known,  since  in  practice  neither  the  signal  polarization  nor  its  correlation 
properties  are  known  a  priori.  Furthermore,  DOA  estimation  algorithms  such  as  [46,  47], 
do  not  make  use  of  prior  information  regarding  signal  parameters.  Thus,  their  performance 
is  likely  to  be  closer  to  that  of  the  CRB  which  assumes  no  knowledge  of  signal  parameters, 
than  to  the  CRB  which  assumes  such  knowledge. 
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We  evaluated  the  CRB  for  some  selected  examples  and  plotted  the  standard  deviation 
of  the  DOA  estimates  as  a  function  of  various  system  parameters.  This  numerical  study 
indicates  that  in  general,  polarization  diversity  improves  the  accuracy  of  direction  finding 
systems  in  the  presence  of  signal  correlation.  In  particular,  we  note  that  the  performance 
of  arrays  which  do  not  make  use  of  polarization  diversity,  is  very  sensitive  to  the  phase  of 
the  signal  correlation,  which  is  a  highly  unpredictable  and  uncontrollable  parameter.  Arrays 
which  use  polarization  diversity  appear  to  be  much  less  sensitive  to  this  parameter. 

In  summary,  the  use  of  polarization  diversity  »  potentially  very  useful  for  direction 
finding  in  the  presence  of  multipath,  provided  that  there  is  a  significant  difference  between 
the  polarizations  of  the  direct  and  secondary  (reflected)  signals. 


2.3  Using  High  Order  Statistics 

Another  area  of  research  undertaken  in  this  project  was  the  development  and  perfor¬ 
mance  analysis  of  parameter  estimation  algorithms  based  on  high-order  statistics,  and  their 
application  to  communications  problems.  Some  basic  theoretical  results  on  this  topic  were 
presented  in  [2,  7,  23). 

High- resolution  direction  finding  techniques  have  been  based  on  the  covariance  matrix 
of  the  received  signals,  i.e.,  on  second  order  statistics.  This  raises  the  intriguing  question 
of  whether  the  higher  order  statistics  of  the  signals  contain  useful  information  about  the 
signal  DOAs,  and  if  so,  bow  could  this  information  be  exploited.  In  [20,  29]  we  made  a 
very  preliminary  attempt  to  address  this  problem  and  developed  a  MUSIC-like  algorithm 
which  exploits  fourth-order  cumulants  to  provide  DOA  estimates.  In  future  work  we  plan 
to  investigate  the  potential  of  high-order  techniques  for  direction  finding  in  the  presence  of 
multipath. 

In  [14,  21,  27,  301  we  presented  some  preliminary  results  on  adaptive  channel  equaliza¬ 
tion  algorithms  based  on  high-order  statistics.  This  work  seems  to  indicate  that  improved 
equalization  performance  may  be  achieved  in  some  situations  by  optimal  utilization  of  the 
information  contained  in  the  high-order  statistics  of  the  signal. 
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